Entangled and disentangled evolution for a single atom in a driven cavity 
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For an atom in an externally driven cavity, we show that special initial states lead to near- 
disentangled atom-field evolution, and superpositions of these can lead to near maximally-entangled 
states. Somewhat counterintutively, we find that (moderate) spontaneous emission in this system 
actually leads to a transient increase in entanglement beyond the steady-state value. We also show 
that a particular field correlation function could be used, in an experimental setting, to track the 
time evolution of this entanglement. 



In a recent, beautiful experiment, Auffeves and 
coworkers Q have verified the prediction Q that, for a 
two-level atom interacting with a single mode of the elec- 
tromagnetic field, in a coherent state with a moderately 
large number of photons, the natural evolution of the sys- 
tem leads to a superposition state in which two different 
states of the atom are correlated with two distinguish- 
able states of the field, in what is often referred to as 
a "Schrodinger cat." Such superpositions of correlated 
states, if pure, are naturally entangled states. Interest- 
ingly, the possibility of preparing such entangled super- 
positions in the above system (which is described by the 
so-called Jaynes-Cummings model, or JCM) arises from 
the existence, in this model, of special trajectories along 
which the joint evolution of field and atom is to a good 
approximation unentangled, i.e., factorizable. It is the 
coherent superposition of such trajectories that results 
in an entangled state. 

The purpose of this paper is to show that a similar 
situation arises in another system of interest, which may 
be easier to prepare and observe in the optical domain, 
namely, a single atom in an externally driven cavity. Un- 
like the JCM, this is an open system, yet, as we shall 
show here, approximately factorizable trajectories exist 
in the absence of spontaneous emission. Moreover, some- 
what surprisingly, we find that the inclusion of sponta- 
neous emission actually helps to create transient entan- 
gled states which are typically more entangled than the 
steady state. We also show that there is a particular field 
correlation function that might be used to keep track, in 
"real time," of the physical processes responsible for the 
evolution of this entangled state. (We note that there is 
a growing interest in exploring the connections between 
"quantum optics" -style correlation functions and entan- 
gled states; see, e.g., Q, and references therein.) 

The system we describe here has been studied before in 
great detail by, among others, Alsing and Carmichael Q 
and Mabuchi and coworkers |5j . The transient regime we 
are interested in here, however, has escaped attention in 
most of these previous studies, because they make a "sec- 
ular approximation" on the Rabi frequency that results 



in an atom-field state that is explicitly disentangled at all 
times (as can be easily seen, for instance, from the ansatz 
in i). An important exception is [fj , where the splitting 
of the field states in phase space, that plays an essential 
role in what follows, was explicitly discussed and illus- 
trated, although the question of entanglement was not 
quantitatively addressed there. (For a very recent dis- 
cussion of how different types of environment monitoring 
of this system may lead to different ways to characterize 
its entanglement, see also 0.) 

The starting point for our analysis is the following mas- 
ter equation for the joint atom-field density operator p: 
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Here, g is the atom-field coupling constant, for the cavity 
field mode, which is described by creation and annihila- 
tion operators oJ and a; tr_ and <r+ are the atom's raising 
and lowering operators; £ is the amplitude of the exter- 
nal, driving field; n is the cavity loss rate, and 7 is the 
spontaneous emission rate. 

As discussed in 8] , in the absence of spontaneous emis- 
sion, approximately unentangled, quasi-pure state trajec- 
tories for this system are obtained whenever the initial 
joint atom- field state is of the form 

|*° ± (r o ,0o)) = ^ (e-^|e) ± \g)) \r e~^) (2) 

where |e) and \g) are the atomic excited and ground 
states, respectively, and |r e _I *°) is a field coherent state 
of arbitrary amplitude r and phase 4>q. Trajectories 
starting from these special states remain approximately 
factorizable and quasi-pure for fairly long times, in spite 
of the dissipation represented by the term k in Eq. Q ; 
they retain approximately the same form as (0, only 
with a time-dependent phase <fi±(t) for the field and the 
atomic dipole, and (in general) a time-dependent ampli- 
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tude r± (i) for the field as well: 

|*±(t;r ^o)) =e i *±W|^(r±(i), < M*))>- (3) 

The overall phase *±(i) will be discussed shortly be- 
low. For 4>±(t) and r±(t), however, we note that con- 
sistency requires that they approximately obey the semi- 
classical equations of motion, derived from by factor- 
ing the expectation values of atom-field operator prod- 
ucts. If one further treats the field as a classical quan- 
tity in these equations, one finds, for £ > g/2 (the so- 
called "strong driving" condition) , a pair of steady states, 
with phases <p u j = T^ss = T sin (<?/2£ ) and amplitude 
r ss = (£/k) cos 4>ss- (The subscripts u and I refer, respec- 
tively, to the "upper" and "lower" steady state, and fol- 
low the notation of The corresponding states of the 
quantum system are respectively 1*°) = |*° (r ss , ~(f> ss )) 
and |*°) = \y°_(r ss ,4> ss )) in the notation of J2J. 

In general, coherent superpositions of states of the 
form (j2J) also need to be considered. Unlike in the or- 
dinary JCM, these superpositions do not remain approx- 
imately pure for as long as the trajectories JSJ them- 
selves, because of the cavity losses; for instance, in Sec- 
tions 3.7-3.11 of it is shown that superpositions of 
|*±( r o,0o)) with the same cf>o decohere at a rate given 
by exp(— g 2 nt 3 /3). Superpositions involving field states 
that start out with very different phases may decohere 
faster, because the photon annihilation operator, acting 
on the corresponding coherent states, will multiply them 
by different phase factors 0. As a result of this, one 
must write down the approximate steady state of the 
system (always neglecting spontaneous emission) as the 
incoherent superposition 

Pss =^° + (r ss , -<l> ss )){^ + (r ss , -0„)| 

+ ||*° (r„, 0a»)>(*i(»-«, (4) 

As an incoherent superposition of product states, the 
state (0J is unentangled, or "separable." 

Consider, however, a single realization of the above 
system, which may have started from a coherent super- 
position of states of the form J5J • Even after the system 
has reached a steady state, and the superposition has 
decohered, it may be argued (at least for as long as the 
individual solutions © remain approximately valid) that 
the decoherence is of the form of a random relative phase 
between the terms of the superposition, a phase that, 
moreover, might be knowable in principle, if we had a 
record of the times at which photons were emitted out of 
the cavity (or alternatively, through a monitoring of the 
transmitted field such as described in [7J). We may then 
ascribe a "conditional" pure state to the system, of the 
form 

|* ss ) = 4=|*° (r„,-0 M )> + ^H*° (r ss ,0 ss )} (5) 



where $' is a random relative phase (time-dependent, in 
general, since the states in the trajectories Q have over- 
all phases that go, for short times, as *±(i) — ^fgrot/2; 
this is analogous to the JCM and is responsible for the 
Rabi oscillations that occur when the two field states 
overlap). For normalization purposes, it has implicitly 
been assumed that the two field states in are orthog- 
onal, which will be approximately the case if r ss sin cj> ss 3> 
1. For large £, this condition becomes g/2n^S> I. 

Using the explicit expressions in shows that 
this is, in general, an entangled state, although not max- 
imally so, since as long as (j> ss ^ the two atomic states 
involved, -^( e *«| e ) + | 5 )), and ^ (e" 4 *- |e> - \g)) , 
are not orthogonal. Assuming the field states are or- 
thogonal, the reduced density operator for the atom 
alone can be written as pa — ije) (e| — | sin cj) ss \e)(g\ + 
%BW</> sa \g){e\ + ±\g){g\, with eigenvalues (I ± sin <f) ss )/2, 
which means that the "entropy of entanglement," E = 
-Tr(pA\og 2 PA) — 1 — cj> 2 ss /2\n2 for small (f> ss . 

Consider now what happens when one has a rela- 
tively small spontaneous emission rate, and a sponta- 
neous emission event occurs. Starting from a state like 
iJSJl, the atom collapses to the ground state \g), so after 
renormalization one has 




= ^+(r ss ,-<p ss )) - |*° (r ss ,-<^ s )>) 

+ ^(|^(W„))H*-(W«)>) (6) 

This expression has been split into two pairs of terms, 
the first one associated with the u steady state of the 
field, and the second one with the I steady state. Each 
one of these is a superposition of the appropriate steady 
state (1*+) in the first parenthesis and |*_) in the second 
one) and another term where the field and atom have the 
"wrong" relative phase. As shown in [g, for short times, 
these nonstationary terms evolve by changing the phase 
of both field and atom at an approximate rate ±g/r ss , 
so the total time-evolved state is of the form |*(t)) = 
^(|*u(*)> +e-**'|*i(t)), with 

|* u (f)) = -^(e- i ^ t / 2 |*° f (r ss ,-0 ss )) 

- e*^"*/ 2 !*" (r ss ,~<f> ss +gt/r ss ))} 

= ^(e-^'^XirssAss -9t/T ss )) 

-e^"*/ 2 |*°_(r ss ,^))). (7) 

Along either one of these two (it or I) branches, the two 
field states involved become approximately orthogonal 
as soon as t ^> 1/g, the JCM's "collapse time." (Note 
that this splitting of the field into four states following a 
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FIG. 1: Solid line: entanglement predicted by Eq. Q. 
Dashed line: the result of a single quantum trajectory calcu- 
lation. Dotted line: the result of a density matrix calculation. 
Model parameters: £ = 0.7g and k = 0.125g. 
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FIG. 2: Atom-field entanglement calculated for a quantum 
trajectory starting from the cavity in the vacuum state and 
atom in the ground state. Model parameters as in Fig. 1, ex- 
cept that the spontaneous emission rate 7 = OAg. The sharp 
(vertical) drops to zero entanglement correspond to sponta- 
neous emission events, that reset the system's wavefunction 
to a state |$)|<;), whatever the field state |$) may happen to 
be at that instant. 



spontaneous emission event was previously discussed and 
illustrated in [||; see in particular Fig. 10.8 there.) If r ss 
is sufficiently large, the phase difference gt/r ss between 
the corresponding atomic states (along the same branch) 
at that time may still be quite small, in which case they 
will still be nearly orthogonal, and the overall state will 
be highly entangled. Specifically, for either state u or I 
we find for small gt/2r ss 

Ecf 1 ( u )-f 2 (u)Bin(2gr sa t)^- - -J- (8) 

r ss 81n2 \r ss ) 

where u = e^ 2 ' 2 / 2 , and fx(u) = (uln((l — u)/(l + u)) — 
ln(l-?i 2 ))/ln4and/ 2 (w) = (2u(l - In 2) + u ln(l - u 2 ) + 
ln((l + u)/(l — it)))/ In 16 are functions associated with 
the overlapping coherent states; at t = 0, fx = and 
ji ~ 0.7, whereas after the collapse time /1 — > 1 and 
f% - > 0, and one has large entanglement provided gt/r as 
is not too large. Eq. (JHJ) applies also to the superposition 
^(|*u(*)> + \^l{t))), regardless of the value of 
as long as the field states in | are orthogonal to 

those in | 

Figure 1 shows a plot of the entropy of entanglement 
E, as a function of time, for the u branch, based on 
the expression J7J for the system's state. The Rabi os- 
cillations actually cause the entanglement to peak some 
time before the collapse is complete. The same figure 
shows also the result of a single quantum trajectory sim- 
ulation (dashed line), and the result of integrating the 
density matrix equations of motion (dotted line) , all for 
the same parameters and initial conditions. The agree- 
ment between the solid and dotted lines indicates that 
the approximation J7J) is indeed quite good. The dashed 
line, on the other hand, suggests that the disruption to 
the relative phase of the terms in JJJ caused by cavity 
losses may sometimes reduce the entanglement obtained 
along an individual quantum trajectory. 

The above analysis shows that, rather surprisingly, 
spontaneous emission may actually help generate sub- 



stantial atom-field entanglement in this system, by peri- 
odically resetting the wavefunction to a state such as iJBJ, 
which can later evolve into something close to a (nearly-) 
maximally entangled state of the form J7J. This expec- 
tation is borne out by further quantum trajectory calcu- 
lations, including spontaneous emission, such as the one 
shown in Fig. 2. Note the pattern: after each sponta- 
neous emission event, the atom- field entanglement nat- 
urally goes down to zero, but then it quickly rises to, 
sometimes, a very high value. This happens even if the 
field state at the time of the spontaneous emission event 
is not the steady state \r ss e ±l ^ ss ) , since a decomposition 
similar to © can be performed regardless of the state of 
the field, and the subsequent evolution will be similar to 
that given in J7J) ; one only needs to use the instantaneous 
phase of the field (rather than the steady state value <p ss ) 
in the \^±) atomic states. 

When the steady state is assumed to be the mixed 
state as opposed to the conditionally pure (JSJ, one 
may assume that the u and I branches exhibited in © 
are superimposed incoherently, although along each of 
them the evolution is still coherent and given by Q . Af- 
ter the "collapse time," when all the field states involved 
have become orthogonal, but the phase shift gt/r ss is 
still small enough to be approximately negligible in the 
atomic states, the resulting mixed state can be writ- 
ten schematically as p = i|* u )(* u | + with 
= -^11) (e'He) - |S» + ||2) (e*«|e) + \g)) and 
|* ; > = i|3) (e-*"|e) + \g)) - §|4> (e^»|e) - | 5 ». (Var- 
ious overall phases have been absorbed in the field states 
|1), |2), |3) and |4).) Treating the field as a 4-dimensional 
system, we find, by the "realignment criterion" [!j that 
this p still describes an entangled state. Specifically, if 
G is the 4 x 16 rearrangement of the 8x8 matrix p, we 
obtain Tr[(GG t ) 1/2 ] = V2 > 1. Thus, we conclude that 
even when the mixed nature of the steady state, for an 
ensemble of identically prepared systems, is considered, 
spontaneous emission does indeed lead to a transient en- 
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FIG. 3: Solid line: h FT (r) — 1) for the same parameters 
as in Fig. 1, based on the expressions 10. Dashed line: the 
result if the pure state I^E'u(t)) (or is replaced by an 

incoherent superposition of |*+) and |*-). 



tangled state, a time of the order of the collapse time 
after the emission event occurs. 

In both the pure and mixed-state cases, the separation 
of the field along each (u and I) branch into a coherent 
superposition of nearly orthogonal states is essential to 
the generation of entanglement. This separation can be 
tracked by using the intensity-field correlation function 
h FT {r) = (a + (0)a g (T)a- (0))/(a + a-)(a e ), where the 
field quadrature operator ag(r) can be calculated, at the 
time r, as ae(r) = U^(T)ag(Q)U(T), and U(t) is the evo- 
lution operator. Experimentally, h FT (r) gives the evo- 
lution of the transmitted field conditioned on the detec- 
tion of a fluorescent photon (i.e., a spontaneous emission 
event) at the time r = 0. The approximation (0 can 
then be used to calculate it; the result (which does not 
depend on whether and \*f?i) are added coherently 
or incoherently, as long as the field states in |\I> U ) are or- 
thogonal to those in \^i}) is plotted in Fig. 3 (solid line), 
where the similarity to the pure-state entanglement curve 
(Fig. (2)) is readily apparent. To better understand this 
similarity, one may consider the following approximate 
result for small gt/r ss : 



h 



FT 



jf 1 / gt 

1 ~ (tan0 ss + usm(2grt))- — 

2r,. s 4 V r.„ 



(9) 



with u = exp(— g t /2), as before. Like Eq. 



Eq. © 

shows Rabi oscillations (although with the opposite sign) 
that die away at the collapse time. In both cases, the 
initial rise of the curves is due to the growing separation, 
in phase space, of the two field states making up the u or I 
branch, although the entanglement eventually saturates, 
around the collapse time, whereas h FT — 1 may continue 
to grow (due to the term tan (j) ss gt/2r ss ) up to a time of 
the order of I /2k. 

We conclude that, through the collapse time, the cor- 
relation function h FT (r) may be used to track the physi- 
cal processes underlying the growth of the atom-field en- 



tanglement in the system, subsequent to a spontaneous 
emission event, although, in order to make the entan- 
glement correspondence quantitative, a fair amount of 
theory needs to be assumed. In particular, note that it is 
not enough to observe an increase in h FT — 1 to conclude 
that entanglement must be growing, since h FT — 1 would 
rise, as a result of the separation of the field states, even if 
the superposition of states making up \^> u ,i(t)} in Eq. (Q) 
was completely incoherent (dashed line in Fig. (3); or 
set u — in Eq. ©), in which case there would be no 
entanglement at all. The Rabi oscillations are thus crit- 
ical evidence that the superposition is coherent and the 
underlying state is entangled. Ironically, these oscilla- 
tions disappear around the collapse time, just when one 
expects entanglement to be largest. In an experimental 
setting, the underlying coherence might be revealed us- 
ing methods such as those suggested in 1], to reverse 
the sign of rotation of the field states and bring back the 
oscillations. 

The feasibility of exploring this entanglement phe- 
nomenon is within reach of current strong-coupled optical 
cavity QED experiments. It will open the strong driving 
parameter space that is different from the one most ex- 
plored to date, the weak driving regime. Further work 
is necessary to expand this to the case of more than one 
atom and to properly account for multiple (partly 
overlapping) spontaneous emission events. 
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